Tema 3 - Dinàmica de fluids 


16 de desembre de 2019 


Que vol dir sa dinàmica de fluids? 

Estudia els fluids en moviment. 

Ara la velocitat si que pot tenir gradients, es a dir s’acceleració local pot ésser diferent de zero i s’acceleració 
convectiva pot ésser també diferent a zero. 

§¥■<>, 


És tot el contrari de tema 2 (Estàtica de fluids). 

En aquest tema estudiarem sa dinàmica en forma integral ( equacions en forma integral o en forma diferencial). 
Nosaltres aplicarem ses lleis de conservació sobretot en forma integral, domes quan sigui necessari les aplicarem 
de forma diferencial. 


Formulació en dinàmica de fluids 


De forma diferencial 
De forma integral ★ 


1 Teorema del transport de Reynolds 


Per definir el teorema del transport de Reynolds abans necessitam fer una serie de definicions i son sobre el nostre 
sistema de control. El fluid és un fluxe de partícules que circulen per un cert conducte adaptant-sa cada vegada amb 
la forma que tengui el conducte, aquest conducte pot ésser des de una canonada, una mànega, un dipòsit, etc.. el que 
sigui. 

Podria ésser que les parets d’aquest conductes siguin deformables. 

El nostre sistema (sistema de control) serà el conjunt de partícules de fluid que en un moment donat ocupen un 
espai, aquest espai és el volum de control. 
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Per tant, per definir el nostre sistema hem de definir quin és el nostre volum de control que un moment donat 
ocupa la nostre massa de control. 



El volum de control ha d’ésser igual o més gran que el volum que ocupa la nostre massa del sistema (. m s i s ). 

Que és el volum de control (V c )? 

Sa regió arbitrària de l’espai que hem triat que en un moment donat és ocupat per sa nostre massa (. m s i s ). 
Important: Sa massa del sistema (m í(í ) és invariant, és una massa que sempre anam seguint. Per tant sa massa 
del sistema és constant. 


m sis = Constant 

Com que hem definit la massa del nostre sistema i el volum de control, podem definir la superfície de control, 
que domes és sa frontera del volum de control, això és la superfície de control. 

1.1 Existeixen tres tipus de propietats. 

• Propietats intensives: Son propietats que no depenen de la quantitat de substància o mida del sistema. Nomes 
depèn de la posició. 
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• Propietats extensives: (N) Son propietats que si depenen de la quantitat de la substància de la mida del 
sistema. 

No te sentit agafar un punt del sistema i dir (i Quin volum hi ha? és una cosa puntual, no te cap sentit. És 
necessari agafar una mida, si agafam una quantitat podem dir que ocupa, sa massa que te serà la seva energia, 
entalpia, moment lineal, etc. 

Anomenam aquestes propietats amb la lletra N. 

N = V, m, u, entalpia,moment lineal, etc. 

• Propietats especifiques: (r?) És tracta de agafar les propietats extensives (N) i dividir-les entre quantitat de 
massa. 

Si N és una quantitat extensiva que depèn de la mida, ho dividim entre la massa del sistema (n? s / s ) i en el 
resultat l’hi direm que és una propietat específica: 

(Propietat especifica) ij = £ = 

- Exemple: Volum específic, Que és? 

És el volum dividit entre la massa o volum dividit entre unitat de mols 

v 5 v 

m N .mols 

Per qualsevol extensiva ( N ) dividit entre la massa del sistema (m s j s ) tindrem una quantitat que ja no depèn de la mida 
del sistema perquè s’ha normalitzat amb la massa del sistema ( m s i s ). 

Sigui N una propietat extensiva i rj la seva específica. 


V = 


N 

m sis 
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Podem obtenir N (extensiva) a partir de rj (específica). 

La propietat N (extensiva) en el nostre sistema serà l’especifica (77) fent l’integral damunt la massa del sistema, 
per definició: 

N s = J, lhis Vdm 

Perquè rj = si cercam N hi ho feim mitjançant integració aleshores: 

N s = fm sls n·dm 

Integram 17 (P. específica) utilitzant la massa del sistema (m í(í ). Però que sabem també, que la densitat és igual a 
un diferencial de massa ( dm ) dividit entre un diferencial de volum (í/v). 

p = ^ si aïllam dm = 4 > dm = p • í/v 

Ara substituïm a l’integral: 

N s = f msis r 1 ·dm 

N s = f Vs . s rj-p-dv (Per definició) 

On ara utilitzam el volum del sistema per fer l’integral. 

• Manera trivial: 
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N S =N 


Donada una quantitat específica (rj) puc trobar la seva extensiva. 

A que és igual el volum? 


(N s )v = f msjs v-dm v (rj) (Específica) 


A l’integral del volum específic v (77) per diferencial de massa (dm) a damunt la massa del sistema 
Després ses integral sa poden fer damunt sa massa ò damunt el volum, és més pràctic fer-ho a damunt el volum 
per això hem fet la definició de volum de control (Vc). 

Nosaltres sabem que densitat: 


O — Uk — massa _ dm . j _ Q j 

r v volum dv ^ am V av 


m 

V 


macro 


dm 

dv 


micro , dv no és el diferencial de volum específic. 


Substituïm el diferencial de massa (dm). 


(Ns)v = fm sis V ■ ^ 

r -1 



dm = p - dv 






dv = que no és diferencial de volum específic 


(N S )V = f Vs v • p ■ dv 




• Exemple: La massa d’un cert sistema (M s ) a que serà igual? 


A l’integral damunt sa massa de (77) i a que és igual 77? Si N és igual a la massa m; 


V = 


N (Estensiva) 

m sis 


Si N = m 


rj = — 


í massa \ 
\ massa/ 


(N s ) M s = f msis 7 ] ■ dm i sabem que 77 = 1 


(N s )M s = f 1 • dm i sabem que dm = p - dv 


(N S )M S = J v p-dv 


(Ns) M s — m s i s 
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1.2 Teorema de transport de Reynolds 


Ara imaginam que tenim un cert volum de control (Vc). 

(Volum de control) 


Sigui sa propietat extensiva (N s ) sa que sigui qualsevol. 

Reynolds el que va demostrar és que: sa variació de la propietat extensiva (la derivada respecta el temps, derivada 
temporal) és igual a dues contribucions. 

• Primera: El que passa a dintre del volum de control. 

• Segona: Si entra ò surt. 

• Exemple: Si tenim massa: sa variació de la massa, si se produeix massa ò sa destrueix (reacció química), si 
entra massa ò surt massa. 

Si formulam el teorema de Reynolds: 



d dt\sis ~ dt ' fffvcV ■ P -dv + ff sc V p -dA 


On 


w\sís = Variació de la propietat extensiva (N)e n el temps, 
j j ■ fff vc 7 ] p dv = Terme local: Evolució en el temps de la propietat del sistema dins el volum de control. 
ffscV P • ^ 'dl = Ei fluxe d’ intercanvi a la superfície de control (fluxe net) 

i d~£ = lt -dA 


El canvi total d’una propietat extensiva ( N ) d’un fluid a l’interior d’un volum de control és igual a dues con¬ 
tribucions; el canvi de la propietat extensiva que pateix momentàniament es a dir la derivada temporal ^ a 
l’interior del volum de control (Vc). 

Si ens fitxam hem de derivar respecta el temps l’integral triple damunt el volum de control (Vc) on: 


/ r) ■ p dv =és la propietat (Ns), propietat del sistema. 
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Les tres integral son perquè és un volum, el volum te tres coordenades (x,y,z). 

Canvi total = canvi momentàniament localment dins el volum de control (Terme local) + El que l’hi esta passant 
a la superfície. El fluxe d’intercanvi a la superfície de control (Terme convectiu). 

Això és el que va dir Reynolds, mes igual el que sigui N; energia, massa, volum, etc.. Varia de tal manera del que 
l’hi passi a dintre si hi ha producció ò destrucció ò a la frontera de l’intercanvi. 

Pot ésser que un dels dos conceptes sigui zero: 


Si domes tenim el segon terme: 


(Terme local) & • fff vc q • p • dv = 0 


(Terme convectiu) ff sc t] • p • • clA ^ 0 


> Problema estacionari 


Hi ha conservació dins del volum de control (VC), domes hem de mirar la superfície de control (SC). 


Ara imaginam que la propietat N és igual a massa per velocitat. 


/V(extensiva) = m ■ V 


Newton va dir que la derivada temporal de massa per velocitat (variació del moment lineal respecta el temps ) 
és igual a la força resultant. Aleshores partint de les equacions de mecànica clàssica amb el teorema de 
transport de Reynolds amb fluid si sa banda de l’esquerra l’hi posam els resultats de la mecànica clàssica i a la 
dreta el transport de Reynolds, podrem trobar força resultant que actua damunt un fluid. 

Coneixent resultats de la mecànica clàssica i 
també termodinàmica clàssica. 

ler Llei de la termodinàmica, conservació de 
l’energ ia._ 

= T.T. Reynolds 

Aquesta mecànica és molt potent. 
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2 Equació de continuïtat 


La variació en el temps d’una propietat extensiva del sistema. 


Afiextensiva) = m =>- rj = 1 (Propietat específica =1) 


í] (específica) = -ff- = ^ (massa) = 1 


Si aplicam transport de Reynolds: 


dN 

dt sis 


TJ ’ p ’ dv ~\~ íí tl·p·V’ ■dÀ 


Terme local 


Terme convectiu 


Í dA = n·dA , 

n —^Sempre és perpendicular a la superfície. 

d A = dx- dy 


En aquest cas aquí trobarem la llei de conservació de l’energia concretament (Equació de continuïtat). 


Que passa si N =massa del sistema. 


Substituïm: 


N — JTlsis = ^’ fi — 1 


\sis — dt ' fffvc fi • P • dv + ff sc fi • p • dA 


dm _ d_ 

dt sis dt 


dt' fffvc 1 ■ p dv + ff sc 1 • p • \f ■ d A 


Aleshores: 


On ^ =0 Per llei de conservació de la massa. 

dt sis 


0 = Jr • fffy C p dv + ffscP-V- d/T Equació de continuïtat de forma integral 


Que passa si N és igual a la massa del sistema {m s j s ) per el vector velocitat ( V"') ? 


N = m sis ■ \f 
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Això és el moment lineal, a partir d’aquí trobarem la conservació del moment lineal. 


• Que passa si N és igual a la massa del sistema per la velocitat el quadrat i partit entre dos? 

N = Wlsis ' ~2 

Això és l’energia cinètica, trobarem l’equació de conservació de l’energia cinètica. 

Això per una banda, d’altre banda alguns resultats general de la mecànica clàssica poden ésser utilitzats per tancar el 
sistema d’equacions, (llei de conservació de masses, segona llei de Newton, etc) 

Per la dreta aplicam el teorema de transport de Reynolds i per l’esquerra aplicam les lleis de conservació de la 
mecànica clàssica així podem tancar el sistema d’equacions. 

• Exemple: La mecànica clàssica diu que la massa ha d’ésser una constant (no sa destrueix ni és pot generar) 

m = cte 


Aleshores sabem que = 0 d’aquí ve ^ \ sis = 0 

•Si N = m s i s ■ V (moment lineal) 

Trobarem que força és igual a massa per acceleració. 


f 


m 


~È 


Però equivalent en els fluids. 

• Que passa si aplicam el primer principi de la termodinàmica? 
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En física general A u = Q + W (variació d’energia interna és igual a calor més treball). Això segons els físics. 
Però nosaltres utilitzarem (enginyers) A u = Q — W . 

És domes un criteri de signes. Si tenim un sistema en el qual realitza treball (W > 0, el treball ha de realitzar 
una tasca) com pot ésser un cotxe, per aquest motiu el calor que entra és positiu i el treball que surt també és 
positiu. 

Continuació: 


•Si N = m s i s que podem extreure de l’equació de continuïtat. 

Aviem dit que si N = m s i s aleshores 77 = 1 i que ^| . = 0 per conservació de masses i així podem tancar 
aquest sistema i el resultat era: 


0 = §i ■ IffvcP · dv + JJscP ■ïf·dï 


{ d~^ = it ■ dA 
dA = dx- dy 


La suma de les dues contribucions és zero. És a dir que si entra o surt matèria, ha d’ésser igual el que s’ha 
produït ò destruït. Això sa coneix com equació de continuïtat de forma integral. 

Que passa si tenim un fluid incompressible? 

Vol dir que la densitat és constant. 


p = cta. 

- Exemple: H 2 O (aigua) però fins a 25 bars. 


Bàsicament els fluids incompressibles son els líquids. En canvi un gas és fàcilment compressible. 

Observació: p = ^ (v =volum especific), la inversa del volum especific és la densitat, i també sabem que: 

P = 7 

Si és un fluid incompressible (p = cta ) aleshores el volum especific tampoc (v) 

Tornant en el fluid incompressible. Que passa si p = cta ? així pot sortir de l’integral i també defora de la 
derivada. 


0 ~ dt ' fffvcP 'dv + IfscP dA 
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Q — P' dt ' fffvcd v + ffsc~^ 'dA 


Sa derivada multiplicada tant el primer terme com el segon, per proporcionalitat podem anular-la (p). 

0 = §~ t ' fffvc^ v ^~ JJsc? ’^A 

On fff vc dv és igual el volum de control (VC). 
fffvc dv = VC 
Convertim l’expressió: 

0 = jj' fffvcdvffsc? dA 

0 = AjV _)_ jj sc f . d Per fluid incompressible 

Per tant zero és igual com varia el nostre volum de control en el temps més l’integral de superfície. 

En aquest moment afegim que el volum de control no sa deforma (parets rígides) no varia amb el temps. 
Podem fer una següent passa: 
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2.1 Equació de continuïtat de forma diferencial: (Si que ens interessa) 

En la conservació de la massa: 

Utilitzant el teorema de Gauss podem convertir integrals de superfície a integrals de volum. Basta afegir l’opera¬ 
dor Nabla ^ V ^ dins l’integral i afegir una integral (integral doble (superfície x,y) passa a una integral triple (Volum 

x,y,z)). 


0 — dt'JffvcP'dv + 


-dl 


sc 


Teorema de Gauss 

(J- 

' IffvcP ' + fffvc ^ ’ (p ’ ’ dv 

Una vegada arribat aquí, ens donam compte que els dos operadors les podem permutar (i p (primera integral)), 
no te res que veure fer una integral de volum a derivar amb el temps. Son operadors que permuten. 


0 = fr 


0 = ///vc[# + ^-(pV) 


•í/v 


a-b = b-a 


Les dues integrals les podem sumar amb una tota sola, si l’integral ha d’ésser zero, l’integrant ha d’ésser igual a 
zero. 


o = M 


vc 


dp 

dt 


+ 


V( P V); 


■ dv 




En forma diferencial 


= 0 


Ara hem trobat l’equació de continuïtat de forma diferencial mitjançant el transport de Reynolds. 

Equació de continuïtat en forma diferencial 


àe+y.fp.yu 


dt 


0 


2.2 Exemple d’aplicació de l’equació de continuïtat en forma integral. 

Hem de ser bastant pràctics a l’hora d’aplicar la teoria explicada anteriorment. 

I anem a resoldre un cas, una simple canonada que hi circula aigua. 

La manera de definir una canonada que hi passa aigua el considerarem com un flux d’aigua que és unidimensional 
(consideram que la canonada te una dimensió) domes té longitud i va passant una certa quantitat de fluid i direm que 
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és estacionari (si feim una fotografia a un temps (t) i després en feim una altre a t +x i veurem que es el mateix). Si 
entren 2 kilograms per segon, han de sortir 2 kilograms per segon indiferentment mirem quan mirem la canonada i 
això vol dir que no depèn del temps. 



Com que és estacionari, el terme local és zero (Jj = 0, les derivades temporals son zero). 
Utilitzant el transport de Reynolds per aconseguir la equació de continuïtat . 

wL = í·fffvc r l·p·dv + ff sc Tl-p-^-dí 


Saben que: 


N — ffï.s/s 
= -AL = Uhií = l 

Wlsis Wl s i s 


*L — dt' fffvc 1 ’ P-dv + ff sc l-p-V-dA 


Com que sabem que no hi ha d’haver variació de la massa (conservació de la massa) el terme de derivada de la 
massa respecta el temps (variació de massa respecta el temps) a d’ésser igual a zero. 


dm 

dt 



Aleshores el sistema quedaria com: 


0 — J;' fffvcP ’ dv + ffscP ' dA 



= lt dA 


Recordam que si és estacionari el terme local és zero. 


0 



+ 


JlscP-f-d^ 


0 
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Per tant queda domes el terme convectiu. 



Ara elegim un volum de control que contengui la nostra canonada, en el nostre cas el nostre volum de control 
tindrà forma cilíndrica. 



Ara hem de resoldre l’integral: 


0 = ff sc p-V-dÀ* \dl = ~rt ■ dA 


A damunt la nostre superfície de control, si aquest cilindre el separam entre la tapa d’entrada, la tapa de sortida i 
el que la part lateral del cilindre i ho representam: 



Ae 





Àrea tub 


Col·loca’m el vector normal de 
cada superfície. 


\y n 



n ▼, n 
n 


? = 


Vector velocitat del flux 




dA= Diferencial normal a l’àrea 
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Hem de fer el producte escalar a damunt cada una de les superfícies 
• Tapa d’entrada: 



Sa component normal (It ) a damunt la superfície. 



Si multiplicam per el diferencial d’àrea. 

V-it ■(-!)■ dA^-V-dA * 


• Part lateral del cilindre: 




I 

n 


V 


\t\ 


It 


cos(90) 



Multiplicar per zero, és igual a zero: 


\?·lt·0·dA = 0 
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0-dA = 0 * 


Tapa de sortida: 



f 

V 


\lt\ 


it 


1 


cos(0) 

1 


Si multiplicam pel diferencial d’àrea. 


V-lt-í-dA^V-dA * 


Aleshores podem aplicar aquests resultats a l’integral: 

0-ffscp.ï.rt 

0 = Jf A ,p·(-V·<iA)+ff A " l p·0·dA + ff A 'P·V·dA 


Superfície de control (SC) 


Superfície de control d’entrada (A e ) 


Superfície de control del tub (A tu / } ) 


Superfície de control de sortida (A s ) 


Si ens fitxam l’àrea del tub (integral A tu } } ) és igual a zero perquè el producte escalar és zero, per tant reordenant 
tot el sistema queda: 


o = SL·,P· v ·dA-SS At p·V·dA 


Assumim que totes les partícules que surten de la tapa van a la mateixa velocitat. 
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Tant la velocitat com la densitat poden sortir defora de l’integral. 

Per la densitat no fa falta considerar que és un fluid incompressible ja que en aquest cas és unidimensional, la 
densitat a tota la tapa té un valor concret, pot ésser així: 


Pe — Ps 

Pe Ps 

La densitat a l’entrada és igual a la sortida 

La densitat a l’entra és diferent a la sortida 


Però damunt la tapa sempre té el mateix valor. 



p e = Densitat d’entrada. 


p s = Densitat de sortida. 

V m = Velocitat mitjana. 

V me = Velocitat mitjana d’entrada. 
V ms = Velocitat mitjana de sortida. 


0 — JXtj Ps ■ Vms ■ dA jj A/ p e ■ V me ■ dA 
0 = Ps ' Vms • ff As dA — p e ■ V me • JJ Ae dA 


ÜA s dA=A s 

ff Ae d A = A e 


^ s 0 — Ps ' Vms ' A s Pe ' Vme ' A e 
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Perquè això sa compleixi (p • V ■ A) ha d’ésser constant. 


Entrada = Sortida 


Pe ' Vme 


' A.e — Ps ' Vms ' A s 


Si ens fitxam allà a on hem posat el volum de control és arbitrari, el conducte si es més llarg s’hagués pogut posar 
a diferent lloc, per tant el producte (p V -A) perquè sa conservi sa massa sempre és igual a constant (sempre). Sa 
llei de sa conservació de la massa a mecànica de fluids s’escriurà d’aquesta manera: 

pe-Ve-Ag = ps-Vg-As 

Per altre banda en el tema 1 varem veure que: sa derivada de la massa respecta el temps és igual a m (cabal 
màssic) [ k g/s] . 


^ = TO i SÍ 

■ 1 

p ■ Q = riï 
[ Q = v m -A 

on Q és cabal volumètric [ m3 /s] 

No ens ha d’embullar, això ens referim globalment no ha un sistema de control. 


Si ara anam substituint a la formula original de conservació de massa a mecànica de fluids tenim que: 


Pe 




Pe ' Qe — Ps ' Qs 




Pe ' Qe — Ps ' Qs 


m P = m ç 


Cabal màssic d’entrada és igual a cabal màssic de sortida. Això vol dir que si entren 4 kilograms per segon (m e ), 
han de sortir també 4 kilograms per segon (m s ). Conservació de la massa. 


• Ara que passa si tenim un sistema on hi ha més d’una entrada i més d’una sortida i totes elles (entrades i 
sortides) les podem considerar unidimensionals. 
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Cada una d’elles (entrades i sortides) les consideram unidimensionals, el fluid domes pot anar per un camí, per 
ell (el fluid) el mon domes te una dimensió malgrat que tot el sistema sigui en 3 dimensions. 

Si consideram tot això i a més que és estacionari, és podria demostrar que el producte (p • -A^j és manté. Sa 
suma de les entrades és igual a sa suma de les sortides. 

E e (p.^.A) e = E ,( p ·^· a) s 




Recordam que la velocitat ens referim concretament a velocitat mitjana sobre les superfícies de sortida i d’en¬ 
trada. Les densitats també son mitjanes. 


En el tema 5 veurem que això realment no és així. 



Si tenim una canonada, sa velocitat del fluid que esta en contacte amb les parets és igual a zero metres per 
segon [m/ s ] per la condició de no lliscament per tant podem suposar que la velocitat tindrà una variació d’un 
mínim a un màxim (V max ). 
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Si feim un tall a la canonada en aquest cas sa velocitat si que depèn de la posició però això és un cas molt 
particular. 

I per això es parla de velocitat mitjana, on: 


V m = (Tema 5) 


3 Equació de conservació del moment lineal. 

Ara direm que la propietat extensiva (N) és el moment lineal. 

N = m sis ■ fi 

On N és igual a la massa del sistema per el vector velocitat j. 
Aleshores sa propietat especifica (rj) és: 


77 = N — m sis' y \ j-j 

m sis m sis 


=fi 


D’altre banda nosaltres sabem que a partir de la segona llei de Newton Çfifi s is = m • ~ct), la derivada temporal de 
la propietat extensiva respecta el temps és la suma de les forces que actuen a damunt el sistema. 

dJLsis _ Y~£ ■ 

dt ~ L r sis 


on 


-fisis = msis ■ V (moment lineal) 

Yfifi sis =m-lï 

~fi sis ^Forces externes que actuen a damunt el sistema, es a dir a damunt el volum de control. 



Llavors per el sumatori de forces externes com que és igual a la derivada temporal també ha d’ésser igual el 
teorema del transport de Reynolds a sa suma de les dues contribucions. 
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Formula de Reynolds: c -§\ sis = §j ■ fff vc r\ ■ p ■ dv + ff^rj • p • • dt 

= t^sis = Jy • fffvcV 'P 'dv + ffscV • p • dA 

d(m s i s -V ) —> 

On <L - d r L· = dl —- i fl (propietat especifica) és V 

lA = = é·fff V ct·p·<lv + fIscV·p·(f·rf) 

Si hens fitxam aquí tenim més d’una equació, son tres, perquè 7^ és un vector. Aleshores s’equació s’ha de 
resoldre per resultant de F x , F y i F- realment son tres equacions, una per cada direcció. Perquè N és un vector, sa 
propietat extensiva te caràcter vectorial. 

3.1 Exemple: 

Ara anam a fer el mateix exemple que abans però en conservació del moment lineal amb un flux unidimensional 
sense saber si és estacionari o no. 



En aquest cas suposant tot el que hem parlat anteriorment del producte escalar. 

El segon terme del teorema del transport de Reynolds: 

Terme convectiu en moment lineal 

Això d’aquí és exactament igual que l’equació de continuïtat però ara ha comparegut un altre terme. 
El que hem dit anteriorment era: 
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JJJoi*'*) 

'-v- / 

Terme convectiu 


ff At f·p·V·dA-Jf Ae V·pVdA 

^·p·V·f{ A! dA-V·p·V·ff Al dA 
V' Pe ■V ms ·ff As dA-\? ' Pe * Vme ' 

‘ Ps ' E mí " Av ‘ Pe ‘ Vme ' A e 


Ara escrivim el teorema de transport de Reynolds: 


dN 

dt sis 


d_ 

dt 



rj ■ p • í/v 


+ 




Terme local + Terme convectiu 


• El terme local: com que no hem dit que era estacionari el deixam com a tal. 


• El terme convectiu en aquest cas és el sumatori del producte vector velocitat 
(V) per àrea (A) a les sortides menys a les entrades. 



per rho (p) per velocitat 


Repetim que el terme convectiu és: 

= Equació de continuitat + el terme ÇV j 

4 L·(pi·VrA i ) s -Z, e (p i ·V r A i ) e 
(és igual a) 

Ara hem de posar a dintre de l’equació de continuïtat la E* on queda: 


(I a (prVrA i ) s -Ze(p i -V r A i ) e )·f 
Ls (V -prVrAi} -^(f -Pí-Ví-A. 


Per moltes entrades i sortides fitxes 


Però el nou terme és el vector que l’hem de multiplicar tant per entrada com per la sortida. Per tant queda per 
un flux unidimensional sa suma de forces externes a damunt el volum de control és igual a aquestes tres contribucions. 

= i-fíí t-p-dv + jjj·p·^·diï) 


vc 


(1 ^ 
sis = J t - JJJ vc ^pdv + l a (f·prV i ·A^ s -l ie (y-Pi-Vi-Ai 


Terme local 
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On si ens fitxam £ ^ V ■ p, • V,- -A,J i recordam que p V A = m (cabal màssic). 

Per aquesta raó podem abrevia les expressions d’entrada i sortida. 

£ Afc = i ■ fffvc f ■ P • ■+ L· ("H ■ c - L· (mr Vi), 
Aquest és el principi bàsic de perquè un coet va per amunt. 



Si estem en tot moment tirant partícules cap abaix, com a reacció perquè s’ha de conservar el moment lineal, si 
hi ha moment lineal cap a baix també ni hi ha d’haver un cap a munt. 

Si ens fitxam, en aquest cas domes tenim una sortida, no tenim cap entrada, per tant surt una força resultant que 
fa enlaira el coet. 


4 Equació de conservació de l’energia 

Una vegada aquí, hem vist l’equació de conservació del moment lineal (dues transparències). Ara donç... 

4.1 Energia cinètica: 

Si Mvariable extensiva) és l’energia cinètica (E c =Variable extensiva). 

N = E c = \·m sis ·V 2 

Per tant la variable especifica (rj = e c ) és: 


4.2 Energia potencial: 

Si N és l’energia potencial (E p =variable extensiva). 
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N = E p = m sis • z-g 


On z [m] és la cota respecta el nivell del mar i g [m/s 2 ] és la gravetat. 
Per tant la variable especifica (77 = e p ) és: 




N 

m sis 


msisjrg 

Wlsis 


= z-g = e p 


4.3 Energia interna: 

Si N és l’energia interna (U =variable extensiva). 

N = U 

Per tant la variable especifica (77 = u ) és: 

11 = ^ 

1 m sis 

Podria tenir altres contribucions, que passa? 

Que la mecànica de fluids ses bàsiques son aquestes. 

4.4 Energia del sistema: 

Continuant, dit això l’energia especifica (r /= <? v /, v ) del sistema és: 

e sis = u + \ - + g ‘ Z 

e s i s =Energia interna + Energia cinètica + Energia potencial 

Aleshores l’energia del nostre sistema (E s j s ) no és més que posar l’especifica (e S j S ) dins de l’integral (extensiva 
en funció de la propietat especifica, teoria anterior). 

Partint que: 

(Extensiva) N s = fff vc fl -p-dv on 77 = e sis 

Substituint: 

Esis = fffy ns (u+\·v 2 + g·z) -p-dv 


u 

m sis 


= = U 
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Una vegada definit això, enunciem el primer principi de la termodinàmica des de el punt de vista dels enginyers 
que és: 

Sa variació de l’energia del nostre sistema es el calor (Q) menys el treball (W). 


E s ís = Q-W [Joules] 


Per què aquest conveni? 

Tot calor que entra és positiu. Per tant si entren en el nostre sistema 2 Joules, el nostre sistema ha augmentat la 
seva energia en 2 Joules. 

En canvi si surten 2 Joules del treball, l’energia del nostre sistema ha minvat 2 Joules, per això es posa un menys, 
perquè l’útil és que les màquines facin un treball. 


w<o 


4 


Q > 0 



Sistema 



Q < 0 



W (Treball) 


w>o 


g(Calor) 


4.5 Conservació de l’energia 


Ara aplicam el teorema del transport de Reynolds. Nosaltres sabem que la derivada temporal de l’energia del sistema 
ha d’ésser igual a derivar la calor ( Q ) respecta el temps menys derivar el treball (W) respecta el temps. 


• • • 


Esis = Q-W 


I aplicant el teorema del transport de Reynolds ha d’ésser igual a: 
T. Transport de Reynolds: 


WL = í • fffvcri p dv + ff sc rj - p \? -dA 


dN I 
dt \sis 


Ès is = Q-W [JoLiles/ s ] 

E sis = Q-W = § J ·fff vc ri·p·dv + JJ sc ri·p·V'·d'í OnJ í l=e sis 


í e sis — M+2'V 2 +g·Z 


Substituïm i reordenam: 


E s is — Q~ W — • fffvcp ■ (u + \ ■ V 2 +g ■ z) ■ dv + ffscP • (u + \ ■ V 2 +g ■ z) ■ -dA j 
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Clar, interessa relacionar calors ( Q ) i treballs (W) amb tot el terme de la dreta. 

A partir d’aquí podem trobar l’equació de Bernoulli, perquè Bernoulli no és res més que la conservació de 
l’energia. 

4.6 Definicions amb treball 

Quan deim treball, hi ha tres tipus de treball dins la mecànica de fluids: 



El treball que interactua amb aquesta massa que hi ha una contribució d’una màquina, llavors una màquina sa que 
sigui que entrega un treball o el fluid l’hi fa un treball a sa màquina. 

Llavors si aquest sistema esta sotmès a una pressió externa, com a mínim la pressió atmosfèrica també realitza 
un treball a damunt sa frontera, també la pot deformar (Pressió = s fff cie ) per tant estarà també la contribució del 
treball com a conseqüència de la pressió que envolta el meu sistema (' W p ). Ens fitxam que tot és treball però no te 
res que veure si l’origen és una màquina o de sa pressió de contorn. I després un altre treball consideram els efectes 
viscosos recordam que en el tema 1; sa viscositat (/i), s’esforç atallant per mu (ji) per derivada de u respecta a 
la derivada de y. 

't — 11 dit 

T — r dy 

Sa viscositat en el final esta fent un esforç i recordam que l’esforç (t) és: 

j7 Forca 

T — _ — -í- 

S Superfície 

Per tant sa viscositat també farà un treball a damunt el sistema, dit d’una altre manera, si la massa del fluid esta 
en contacte amb un altre fluid, existeix un fregament pel fet que hi ha una viscositat, aquest fregament farà un treball 
a damunt el sistema. 
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Aleshores quan parlam de mecànica de fluids, hi ha tres tipus de treballs: 

Màquines (W m ) 

W < Pressió del voltant (W p ) 

El fluid en si mateix per viscositat (W p ) 

Donat que el transport del teorema de Reynolds ens sortien derivades temporals, aleshores posam punts per tot. 

w = w m +w p +w l·l 

I direm que és el treball dividit entre el temps realitzat en el volum de control (VC) per una màquina on hi ha 
efectes de treball, per pressió i de viscositat. 

Això és una simple definició, ens interessa desglossar el treball amb aquest tres termes per poder arribar a l’equa¬ 
ció de Bernoulli (un dels motius). 

Si hem separat el terme de treball en tres contribucions, cada una la podem escriure ò definir de la següent manera: 

4.6.1 Treball per pressió: 

W p = ff sc p(^·‘l^) 

Anam a deduir-ho: 

Nosaltres sabem que el diferencial de treball que esta actuant a damunt una superfície com a conseqüència d’una 
pressió, ho podem escriure com sa força que fa la pressió per el diferencial de / (dl). 

dW p = d/È p ■d l 
Treball = Forca ■ Distància 

El producte escalar és necessari per si el sentit (direcció de F) és perpendicular al diferencial de /, aleshores no 
sa fa treball. Però clar nosaltres volem la variació temporal (dW p ), llavors derivam respecta el temps i el diferencial 
de 1 (d~t) sa converteix amb un diferencial de velocitat. 



Aleshores: 
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dw p = df t 


dW p = df p 


’Wn = d~f n ■ df 


D’altre banda; que és 7^, 


Sa força de pressió no és més que pressió per diferencial d’àrea. 


~f p =P-dÍ 


Si ajuntam les dues parts: 


dW p = dfp ■ df 
f p = P·dt 


dW p = d~f p ■ df 


dW p = P ■ d~f ■ df 


I després integram a damunt la superfície: 


w p = 


SSscP- (?• 


Ara recordam el primer principi de Termodinàmica: 


Calor — Treball 


E sis = Q-w 


Però el treball l’hem separat en tres termes. 


• • / • • • 


E sis = Q-[W m + W p + W l· 


E*is = Q — W m — Wp — Wp 


I el terme de pressió l’hem posat en forma d’una integral. 


Wp = SJsc p ■ 


Aleshores: 
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Es is — Q-W m -fJ sc P·(ï?·dí)-W fl = A que ha d’ésser igual? Al transport de Reynolds. 
E sis = Q-W m -ff sc P·(f·d^ y )-W tl = Reynolds 
E sis = Q-W m -ff sc P·(y·d~l)-W l·l = §- t ·fff vc r 1 ·p·dv + ff sc r,·p·^·d^ 


I sabíem que sa propietat especifica era: 


V — Csis — u + \ ■ V 2 + g ■: 


E] is = Q-W m -ff sc P·(V·d~í)-W l·l = £ t ·fff vc p·(u + \·V 2 +g·z)·dv + ff sc p·(u + l·V 2 +g·z)·(V·d'l 
Si ens fitxam, el terme de pressió de l’esquerra (l’integral) pot passar a l’altre banda sumant. 

E] is = Q-W m -W l·l = j i ·fff vc p·(u+ l ï·V 2 + g·z)·dv+ff sc p·(u+ í ï·V 2 + g·z)·(y·d~l)+ff sc P·(V·d'2 


On 


JfscP ' (u+\-V 2 + g-z) ■ (J-dA) +ff sc P· Çï^-dA 


Tenen el mateix integrant i sa mateixa integral (propietats d’integrals) 


E' sis = Q-W m -% = j r fff vc p·(u + ±·V 2 +g·z)·dv + ff sc p·(u + ±·V 2 +g·z)·(v'·dl)+P 


Si feim el tros de 


p·(u + í·V 2 +g·z)+P 


p·(u + í·V 2 +g·z)+P 

4 

p -u + p - ^-V 2 +p -g-z + P 

4 

p·u + p·j·V 2 + p·g·z+j 

Ara hem de treure factor comú, però ara ho farem d’una manera molt especial, ja que volem que P tengui també 
p, aleshores 

1 t/2 


Ara treim factor comú: 


p·u + p·±·V 2 +p·g·z + p·^ 


p ■ [u+^·V 2 + g·z+^ 


Una vegada arreglat aquest tros d’expressió podem posar l’expressió completa: 
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E sis = Q-W m -W l·l = j- t ·fff vc p·(u+\·V 2 + g·z)·dv + ff sc p·(u+\·V 2 + g·z)·(V·d^+P 

A, = 6-W m -W At = |·/// vc p·( M + l·7 2 + g·z)·Jv + // 5C p·( M + g + l·7 2 + g·z)·(y·jl) 

Això ja és un resultat més general, per què ens interessa ficar la pressió? Perquè la pressió és una propietat del 
fluid i ens interessa. 

Sa pressió que actua a damunt el sistema vol dir que el sistema te la mateixa pressió (pressió que esta a damunt 
sa frontera). 

• Observacions: 

u + ^ és entalpia especifica. 

La inversa de la densitat és el volum especific. 
h = u + ^r = u + P v = h 

r 

H = v + P-V 

h = u + P • v (Especific) Pv = £ 

r 

Continuació: 

Ara imaginam que tenim un flux estacionari; 

Sa part local és zero perquè qualsevol derivada temporal serà zero. Per tant queda la següent expressió: 

Q — W m — W j _ l = j l · JJJyc p ■ (u + \ -V 2 + g ■ z) ■ dv + ff sc p ■ (u + ^ + \ -V 2 + g ■ z'j ■ -dA^j 
§1 ' IffvcP • (w + \ ■ V 2 + 8 ■ z) ■ dv = 0 (Descartam la part local) 

Per tant per fluxe estacionari l’expressió és: 

Ò-Wl-^=//5Cp·(^+g + l·V 2 +g·z)·(^·^) 
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Si és estacionari i tenim entrades i sortides fitxes unidimensionals. Aplicarem tot el que hem dit abans (producte 
escalar). 

En el final nomes ens interessa aquestes variables, que és el que els hi passa en els valors que prenen a la sortida 
com ha valors mitjans i també a l’entrada. 

Q-W m -W l·í =J!scP·^ + ^ + l · v2 + 8 ·z)·(^·d^ 

■D- [m = p·Q] 

Ò-W w -l^=i: í m·(u + ^ + i·y 2 + g·z)-L e m·(u + ^ + i·y 2 + g·z) 

Quan és unidimensional nomes ens fitxam el que passa a les sortides i a les entrades. Això comença a semblar-se 
a l’equació de Bemoulli. 



Tenim l’equació d’abans però sense sumatoris, on la sortida és l’index 2 i l’entrada l’index 1. 
El que nosaltres coneixem es la velocitat mitjana, l’àrea, la densitat i la cota. 

Q-W m -W l·l =rh 2 ·(u2 + f 2 +\·VÏ + g·Z2)-m l ·{u l + ^ + \·Vl+g·zi) 

Aquí ja començam a veure coses de l’equació de Bemoulli. 

Que passa ara? Apliquem la conservació de la quantitat de massa. 

m.2 = ríi\ = th 
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El cabal màssic en el punt 2 ha d’ésser igual que el cabal màssic en el punt 1. 


Aleshores si teníem que l’expressió era: 

Q-W m -W tl =m 2 ·(u2 + ^ + ^·V} + g·z 2 )-m 1 ·(u 1 + jl + í·V 1 2 + g·zi) 

Si els cabals màssics son iguals ( m 2 = rh\= m ), les podem passar l’altre banda sense cap problema dividint. 


Q-W m -W l·l =m2·(u2 + % + \·VÏ + g·z 2 )-m 1 ·(u 1 + ^+ 1 1 ·V?+g·zi) 
Q-W m -W l·i =m·(u 2 + ^ 2 +\·VÏ + g·z 2 )-m·(u l + l· { + \·VÏ + g·zi) 



- Resum de classe anterior: 

Aplicació de teorema transport de Reynolds a la conservació de l’energia. 

Definir el criteri de signes on el calor absorbit és positiu, el calor que sa cedia cap a fora del sistema és 
negatiu. El treball que sa realitza a fora és positiu, i el treball que s’inverteix cap a dins el sistema és 
negatiu. 



Sa variació d’energia respecta el temps és igual a sa suma de dues integrals pel transport de Reynolds 
El treball el separàvem en tres termes: 


* Treball mecànic (motor). 

* Treball de tipus viscos. 

* Treball per pressió que envolta el fluid. 
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Exemple una bomba centrífuga o una turbina eòlica, aquest motor farà un treball a damunt el nostre 
sistema, i direm que és un treball mecànic. 



El mateix sistema te el fluid que te viscositat, aleshores te un fregament d’ell mateix i amb un altre fluid 
que pugui estar en contacte, aleshores el treball del fluid 1 a damunt el fluid 2 seria de tipus viscos. 
Després a part, el sistema pot estar envoltat d’una certa pressió que el pot empènyer cap una banda o cap 
a l’altre seria el treball. 

Recordam que: 

W = F S (S = espai (desplaçament)) 

Aleshores el treball interessa separar-ho amb aquest tres fenòmens, perquè son absolutament diferents. 


W = F ■ l 



4 4 

dW _ F çU 

dt dt 
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[J/s = watt] W,„ = Treball derivat respecta el temps = potència. 

Conservació de sa massa m\ = m 2 

Cada vegada que agafam absolutes extensives i les dividim entre quantitat de massa ens queda valors 
específics 

u = m (« = Energia interna) 


Calor i treball intercanviats per sistema per unitat de cabal màssic. 


Q-W m -Wa 


Q=a 

m 


W m 


= w, 


w„ 


rm 


= W L 


q — — — (ü2 + ^ + \ • V r 2 2 +g·Z2^) — («1 + ^ + \ ' V\ + <? 'Zl) 

Partint de l’equació anterior, passam a l’esquerra la informació del punt 1 + \ • V\ + g • zi) menys el terme 

d’energia interna ( u \). 

+ \ • V? + g • Zi + q — w m — w ll = (ü2 + ^ + \ • E 2 2 + g • Z2 ^ — mi 
P assam la part de potència calorífica ( q ) a la dreta i deixam la informació del punt 2 («2 + 0 + \ • E 2 2 + g • Z 2 )- 

+ 5 • Ej 2 + g • Zi — H’m — ^ + 2 ' ^2 + £ ■ Z2 ~ <? + «2 — «1 

Igualment passam a la dreta el terme de treball de viscositat (w^) i el treball mecànic (w m ) queda igual. 

^ + \ • V\ + g • Zl — W OT = ^ + 2 • E 2 2 + g • Z2 — q + M 2 — M 1 + w l·i 

Ara si ens fitxam el conjunt —q+ «2 — «i + wy, (calor intercanviat quan el fluid viatge per la canonada + variació 
d’energia interna del punt 1 en el punt 2 + el treball de viscositat, un fregament entre ell mateix) 

és igual ag-l·ii. 


-q + u 2 -ui+w l· i = g-h L 


Tot això se’l coneix com a pèrdues per fregament entre les parets i efectes viscosos (Pèrdues per fregament). 
Tot això s’agrupa i és un sol terme {ghi), per què? 

Perquè ficar-nos a cada cosa és la teoria de diversos temes. Però com que sa pot agrupar i mesurar amb un 
experiment i tenir un sol valor, ho feim així i és molt més pràctic. 

Aleshores substituïm aquest nou concepte i l’equació de Bemoulli és: 
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^ + \ • + 8 • z\ — w m — ^ + \ • + § • Z 2 + g • 

Aquí hi ha dos termes nous (w m i g • hjJ. 

Energia punt 1 = Energia punt 2 


Però que passa si entre el punt 1 i el punt 2 hi tenim una màquina? 

Com que és un motor, perquè el motor el que fa és agafar aigua del punt 1 hi ho puja en el punt 2, llavors aquest 
terme és w m , com que anam del punt 1 en el punt 2 i no és un fluid ideal (te viscositat) existeix un fregament 
entre ell mateix i amb les parets. 

És a dir l’energia que tenim en el principi més l’energia que l’hi dona una màquina és l’energia que te en el 
final més unes pèrdues (les pèrdues poden ésser negatives o positives). 

Les màquines que interactuen damunt el sistema poden ésser de dos tipus: 

- Donen treball en el sistema. 

- Absorbeixen treball. 


Coses a tenim en compte: 


W m = m • w m < 


W m 


■É 


= m 


[J/s = watts] = 


w„ 


[kg/s] 


■É 

[J-kg] 


m ■ w m = W m 


4 


kg . J_ _ 'k&·J _ l 

s kg s-kg s 


= 1 watt 


W m 


Variable 

Unitats 

Definició 


J 

watt \ 

kg _ kg_ _ J_ 

S X kg 

Treball mecànic 

• 

m 

kg 

s 

Cabal màssic 

W m 

j = 1 watt 

Potència 
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• • 

Si W m és negatiu (W m < 0) es tractarà d’una bomba centrífuga (bomba). Capítol 6. El treball és negatiu perquè 
és un treball que l’hi dona a l’aigua. 

Si és una turbina eòlica per exemple, el terme W m serà positiu perquè extreu energia del fluid. Surt energia en 
forma de treball a partir del fluid (W m > 0). 



• ^ 

Si w m és positiu W m també ho serà. I m és un cabal màssic [ k s/s]. 


Imaginam que tenim un fluid que transportam d’un punt 1 a un punt 2 i posam un termòmetre tant en el punt 1 
com en el punt 2 i encenem una bomba que impulsa aigua. La temperatura del punt 1 serà igual que en el punt 
2 ? 


No, per un terme de viscositat (fregament) podria augmentar la temperatura. L’increment de temperatura pel 
fregament en si mateix serà menyspreable però si que sa dona. La bomba també cedeix un calor en el fluid. Un 
altre terme és l’energia interna, aquest sobretot depèn de la temperatura, és a dir que si aportam calor, variarà 
la seva temperatura. 


Com veim és un sistema bastant complexa analitzar exactament que l’hi passa a l’aigua, però en aquesta 
assignatura de màquines hidràuliques la hipòtesi és que sa bomba no aporta calor o és menyspreable. Si no 
aporta calor, la seva temperatura no canvia i per tant la seva energia interna tampoc. 


q 


Aquests tres termes seran pràcticament zero < M2 


U\ 


El que no serà zero serà w ^ perquè els fluids no son ideals, tenen viscositat. 


Per què és diu g-hil 

Perquè en el final és una definició: g-hL = 4 =>• h^ = \ Que és equivalent com a pèrdues d’energia potencial. 

O 

hi És una alçada. 
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Direm que: ghi = 4 el que sigui, és una definició. Que vol dir? 

Que si agafam aquest quatre i el dividim entre g, tenim una alçada (hi) equivalent com a perdeu d’energia 
potencial. 

g ■ h L = 4 => h L = ~ s [i m.c.a .] 

És per assimilar sa perdeu per fregament, el que en el final fa que no estigui en 7,2 sinó que estigui un com més 
amunt perquè hi ha hagut aquesta perdeu. 

Dit d’una altre manera, per arribar a zi sa energia del punt 1 ha d’invertir amb aquestes pèrdues per arribar en 
el punt 2. 

Així és més pràctic. 

Quasi sempre el terme de màquina (W m ) en el tema 3 serà zero, perquè no hi ha cap motor. 

Sa perdeu es pot expressar com: 


Pl = P ■ g-hL [Pascals (Pa)] 


PL _ P-g^£ 
P P 


= g-h L 


Si ses pèrdues de fregament son, exemple 4 Pascals, vol dir que el terme ghi = p ' 8 n hL (s’ha de dividir entre p) 
Una pressió sempre sa pot escriure com a p ■ g-hi (pressió estàtica) que la seva unitat son Pascals [Pa]. 

Llavors aquest és el terme que introduïm a l’equació de Bernoulli però hem de pensar que l’equació de Ber- 
noulli nosaltres la podem plantejar “així com nosaltres colguem" 


[gh L \ = Pa /kgm 3 =^j 

- Exemple: 

l·iL = 2 m = 2 [m.c.a.] 

Si hi = 2 m aquests dos metres, son dos metres de columna d’aigua [m.c.a]. Aleshores que val ghiP. 
h L = Im.c.a. => gh L = 9.81 • 2 

En un exercici s’ens pot demanar trobar ses pèrdues de pressió a conseqüència del fregament. És tracta de 
trobar ghL i multiplicar-ho per la densitat (p) 


gh L -p = [Pa] 
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4.7 Equació de conservació del moment lineal en forma diferencial i el tensor d’esforços. 

Un esforç no és més que força partit per superfície. 

—y 

x = j és molt similar a la pressió: P = j 


Quina diferencia hi a entre un esforç i una pressió? 

La direcció que nosaltres feim la força. Sa força és realment un vector. Per tant en funció de quina direcció te, 
posam un subíndex. 

En el cas d’un fluid newtonià varem veure que l’esforç que fèiem era a damunt de la cara normal de l’eix y en 
direcció x. 



El primer subíndex vol dir que és el vector normal. 




Segon Subíndex 



El segon subíndex vol dir la direcció de l’esforç, en aquest cas a damunt l’eix x. Aquest esforç se l’anomena 
ateant. 

Altres esforços poden ésser: 
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A tracció 






A compressió 


A torsió 


Recordam que en el capítol 2 l’equació general d’estàtica de fluids és podia expressar com: 


4=-^p+pt 


dv = Diferencial de volum 


Però aquí en aquest fluid si recordam del capítol 2, quin son els efectes que varem posar? 



Teníem un cub en el qual hi havia les forces de pressió i la força de gravetat (capítol 2). El cub no sa movia. Que 
passa quan el fluid esta en moviment? 
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Quan es mou apareix un esforç ateant, per tant ara a part de l’estàtica que actua sa pressió i sa gravetat aquí hi 
haurem d’afegir un altre terme en el cas que el fluid estigui en moviment. 

Quan es mou hi ha una inèrcia que s’oposa en el moviment. 

En els sòlids rígids és la massa (massa inercial) 

En els fluids és la viscositat. Com s’oposa? 

De la manera com va dir Newton. 

V = M ‘ % (capítol 1) 

Aleshores aquest és el terme que hem d’introduir aquí: 

Ara ja no és tan senzill i anam-ho a justificar com a partir de que un fluid ja no és ideal ni és estàtic, resulta que 
és viscos i afegim el terme anterior. Dit d’un altre manera la força que aplicam a damunt el fluid ha de contrarestar 
les forces viscoses perquè s’oposen en el moviment. 

Sa recorda (capítol 1) que la definició de fluid Newtonià: Que teníem a l’eix y, per tant aquest esforç actua a 
damunt sa cara que esta a la normal de y però en direcció x. I aquest esforç és igual a ”mu” (/i) per derivada de u 
respecta a y. 



On u és la velocitat en la que sa mou la superfície. 

Però ara que passa? 

En el nostre cub que nomes tenia forces de pressió i de gravetat. 
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Però ara poden actuar a damunt de qualsevol cara i fer-li un esforç que nosaltres vulguem. 



En el cas d’una torsió seria una suma. 


Recordam que esforç és ^ 


Esforç en direcció j sobre la cara normal a l’eix i. 


Si és tracció serà Zu 


Si recordam en el capítol 2 que teníem una pressió en el punt x i a l’altre banda x + dx, en el final donava igual si 
una banda teníem 100 i a l’altre banda 100 + 1, el que és important és la diferencia (|j). 






I 

I 

I 

: 

i 

i Px+5Px 

: -v 

i 
i 
• 
a 


dx 


X 

100 [Pa] 


x+dx 

1004® [P I 



És la diferència 


El que afecta és la derivada. 


P-(P + SP X ) = SP X 


Que s’hauria de fer per saber l’esforç resultant en direcció v del dibuix del cub d’esforços? 


y(j)‘ 

z(k)>/ 


res yx j xx i zxk 
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Suma de totes les components de direcció x. Agafar el fluid i du el fluid en direcció x. 
Una vegada que tenim això, és molt pràctic calcular l’esforç resultant de cada direcció. 
Així podem definir el tensor d’esforços com: 


Tij = 




^XX 

^yx 

^zx 

^xy 

Tyy 

T ’zy 

T'XZ 

ïyz 

% 


=> en direncció x 
=> en direncció y 
=> en direncció z 


Recordam del capítol 2 que el que realment importa son el gradient dels esforços ( 8P X = — yjy) que son els que 
originen una força neta sobre sa superfície de control. 



Domes tenim en compte els afectes de pressió, que ens va sortir? 

Que el diferencial de força per unitat de volum és igual a menys el gradient de la pressió (recordatori tema 2). 


dFx = P ■ dy ■ dz — (P + dPx) ■dy-dz 


P + 8P X (8P x = dPx=^-dx 




P+f x dx 


4 


dFx = P ■ dy ■ dz — [P+^f-·dx] ■dy-dz 


4 


dFx = P -dy-dz— [P ■ dy • dz~\~ ^ • dx • dy ■ dz 
dFx = — ^ • dx-dy-dz 


dFx 

dx-dy-dz 


dP 

dx 


4 


dFx _ _dP 
dv dx 
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diferencial de foga 
diferencial de volum 


=>- on és el gradient de pressió (Capítol 2) 


Si això ho tenim en ment podem fer exactament el mateix però si abans teníem una sola component en aquest cas 
en tenim tres. 


dP 

P + • dx ■ dy ■ dz 

dx 


4 


(Capítol 2) 


^x + ^-dy^-dx-dz 

Aleshores és fer la derivada z xx respecta a x i les altres: 


=4> • dx-dy-dz 

=>• ■dx-dy-dz 

=>- •dx-dy-dz 

I sumar-les, abans domes hi havia l’efecte de la pressió, però aquí hi ha tres elements que fan el mateix fenomen 
que du el volum cap a la dreta. 


• Repetició Cap. 2 (Amb pressió P en direcció x) 
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Amb un sistemes de coordenades cartesianes unificam els 3: 

d ^ = -(% + % + %) • dx-dydz 

^ = + + -dx-dy^dz 

• Feim exactament el mateix amb l’esforç (t) 


En pressió 

(. P x ) P-dy-dz- (P+jy- dxj ■dy-dz 


Exporta el mateix però amb esforços 


Si nosaltres suposam així com abans suposàvem que hi havia una funció ( P ) que era coneguda, aquí suposam 
que tenim una funció d’esforç (t) desconeguda. Hem de treure l’esforç net (£t) per cada cara del cub, en 
aquest cas hem tret l’esforç net si l’esforç és fa en direcció x. 

- Cares amb component normal en el vector y: 

Aleshores deim (per un criteri de signes) que la fletxa de sota deim que va cap a l’esquerra i la de dalt cap 
a la dreta. 
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(jyx + ■ dy^j ■dx·dz+ (~Ty X ■ dx ■ dz ) 

+ lïf ' dy-dx ■ dz-'Vy^xL·L·· cU 


=> -gy ■dx-dy-dz 


- Cares amb component normal en el vector x: 



(jxx + • dx) ■ dy ■ dz + (~r xx ■ dy ■ dz) 

'r^rdy^dz + yjf-dx-dy- dz~^ryr~thu4z 


=> yjy ■dx-dy-dz 

- Cares amb component normal en el vector z: 
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( Z zx + ^f ’ dz ) • dx ■ dy + (-Tzx ■ dx ■ dy) 

^P r d^dyy\- yyr dz-dx■ dy^X^r-dx^jiy. 

=y yjy ■dx-dy-dz 

- On recordam del capítol 2: 

On dv és un diferencial de volum. 

Forces viscoses T yx = d'% 

On T yx és l’esforç en la cara normal y i en direcció x. 

L’esforç aplicat a cada una de les cares pot tenir tres direccions diferents. 

La resultant de l’esforç en direcció x és (vector). 

-4 —ï. 

t — txx ' 1 "L tyx ' ] "F tzx ' k 

dF x = — • dx- dy ■ dz+ yjf- ■ dx- dy ■ dz + yyy - dx ■ dy ■ dz + yjy ■dx-dy-dz 

On — ^ • dx ■ dy ■ dz ho sabem del capítol 2 i ■ dx-dy ■ dz + ~yy ■ dx-dy ■ dz + yyr • dx-dy-dz és el que ens 
ha sortit en el resultat per de l’esforç de cada cara del cub en direcció x. 

Sabent que dx ■ dy ■ dz = dv és el diferencial de volum i ho podem passar a l’altre banda dividint. 
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dF x _ _ dP i dx xx i dTyx | dz Z x 

dx-dy-dz ' dy ' dz 


dF x __ d_P <9t xjc dXyx , <9t zjc 

dv dx ' dx ' dy ' dz 


Si recordam en el capítol 2, i que si teníem una partícula de fluid representada com un cub en el qual el seu 
entorn nosaltres coneixem sa pressió que era una funció de l’espai i a més de temps. 


P = f(x,y,z,t) 


Però com que en el capítol 2 era estacionari, no depenia del temps. 


P = f(x,y,zÀ)=>P = f(x,y,z) 


Que vol dir? 

Que nosaltres coneixem la pressió que és una funció i ha d’ésser continua i derivable per l’hipòtesi del continu. 
Coneixem el valor de la pressió. 

- Exemple: En el punt v la pressió serà P x , si nosaltres ens movem en un diferencial de x. 



40 — 40 + x = x (Alguna cosa) 

Domes queda el terme de la derivada. 

En el terme de l’esforç, és que no domes actua a damunt una sola cara, si no que la component x tenim que pot 
actua a damunt tres cares. 

Però un cub te sis cares aleshores tindrem un conjunt de cares. 
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Laterals Alt i baix Davant i darrera 


Així com aquí el criteri de signes és que la pressió va encontra la superfície, el criteri de signes pels esforços és 
que si la cara de baix va cap a l’esquerra, a d’alt deim que va cap a la dreta. Però domes és un criteri de signes. 
Una cara amb l’altre sa cancel·leran i domes afectaran les derivades. 

Esforç (t) = j 

On dx-dy = ds és un diferencial de superfície. 

Senzillament és un balanç de forces que estan actuant. 

Com que nomes sobreviu el terme de les derivades, podem escriure que la força resultant en direcció x serà: Les 
tres noves components que han sortit respecta el terme de la pressió que ja teníem del capítol 2. 

(^ — 'j^·dx· dy ■ dzj i ara afegim les tres direccions (tres components) que estiren el fluid cap a la direcció x com 
a conseqüència d’un esforç que esta actuant a damunt el fluid. És un altre fenomen. 


Fenòmens 

Efectes de pressió 


Efectes d’esforços que estan donats per la viscositat pròpia del fluid que l’envolta (ell mateix). 

El fluid interactua en si mateix i amb altres. 


Si ara feim el mateix que en el capítol 2, ens sortia un diferencial de volum ( dv ) que és comú a tots els termes. 

dF x = — ^ - dx-dy-dz + - - dx■ dy ■ dz+ ■ dx-dy ■ dz + - dx-dy ■ dz 

Passam el diferencial de volum ( dx -dy-dz — dv) dividint. 

4Ilx _ _ d]L I FZxX I d'Tyx | à T- y 

dv dx'dx'dy' dz 

Que després afegim el terme de la gravetat (no surt ara). 

Ara afegim el terme dels esforços que si ens fitxam: 
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d T X y | dZyx . d X zx 
dx ' dy ' dz 


És una operador Nabla (V) a damunt el vector en direcció x que seria la primera filera del tensor d’esforços. 
Aleshores podem abreviar l’expressió: 


dF x _ _ dP | dz xx i d T yx | dz zx 
dv dx'dx'dy' dz 


dZjçx i dlyx dz zx 
dx ' dy ' dz 



Per tant: # = -f+V.^ 


Això nomes és en direcció x. 


Ara hem de fer exactament el mateix però en totes direccions. 


• Tenim el balanç de forces en direcció x: 



- En direcció y, que si ens fitxam hi ha el terme de la pressió, el terme de esforç de la viscositat però a més 
en direcció y també hi ha la gravetat. 

(eix)t y ^ = g-(-T) 

* D’on surt ^ = — jy + p ■ g+ V • ? Anam a resoldre-ho: 

Suma d’esforços a l’eix y: 


T—= T + T -> + T 
''ress miT ' ™ / ' % 


yy j 


xy i 


zy k 


• T -y: 

yy j 


• T 

xy i 


zy k 


(%• + ^ ^ ■dx·dz+ (-% • í/x• Jz) 
^ • dy-dx ■ dz =r Zy X -dx^dz L· 


z -> = ■ dx-dy-dz 

yy j ày j 


(zxy + • c/xj • dz ■ dy + (~z xy ■ dz ■ dy) 

d T 

+ ~df dx dy■ dz^Zyy-dz^dx. 
T xy -? = d -£-dx-dy-dz 
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Sumam les components y : 


(% + ' dzj • dx ■ dy + (—% • dx ■ dy ) 

d'ï 

^c^dx^dyF -yy ■ dz-dx- dy^T^y-dxzzdy 
T -> = • dx-dy-dz 

zyk dz ■ 

^ yy; i zy k 

^f-dx-dy-dz+^-dx-dy-dz+^· dx-dy-dz 

Ara hem d’afegir un terme del tema 2: ^ = — V •P + p·"^’ 

Aleshores afegirem : — + p • 

El diferencial de la força en direcció y és: 

dF y = +P ■ it + W 

dF y = +p -~^ + -dx-dy-dz+^f- -dx-dy-dz+^f -dx-dy-dz 
dFy = -^+p-~Í + + ^ + -dx-dy-dz 

dF y = — + P • ^ ^ ' T xy + ^ ' T z;y) • dx-dy-dz 

dF y = - + p • ~f + 7 • -dx-dy-dz 

dx-dy-dz ~ dy + d 8 v J 

dv ~ dy + P 8 + v T y * 
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Torsió 
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En direcció z: 


Si feim el mateix per treure el balanç de forces en direcció Z: 



y 

E 

c 

( 

+ 

Z _ /• / 


3i A 

1 



!,/ 


r 




- - -i 


^ . :T.:: 


Suma d’esforços T z 
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* Suma d’esforços a l’eix z: 


T Tà z : 


T t> + T T -7» 


zz k 


xz i 


yz k 


• T —y ! 

xz Z 

(jxz + ^ • dz ■ dy+ (~T XZ ■ dz ■ dy ) 
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■dx dz ■ dy^X^-dz^bt 


x -f = ^jf ■dx-dy-dz 


yzj 


T yz + —jyy • dy^j -dx·dz+(-T yz ·dx·dz ) 

• dy-dx- dz-^C^y-dxzdZy 


T -> = -dx-dy-dz 

yz i ay J 


zz k 


x zz + ■ dz^j -dx-dy + (— T zz • dx■ dy) 

^z^dx^dXjV ^jf -dz-dx■ dy =r x^dxzzdy_ 

T t> = ■dx-dy-dz 

zz k dz 7 


Sumam les components z 


W ^ + T * Z ? + T zz f 

T^z = ^ -dx-dy-dz + -dx-dy-dz+^f -dx-dy-dz 

Ara hem d’afegir un terme del tema 2: — ^-dx-dy-dz 
El diferencial de la força en direcció z és: 


dF z 


dF z — + ^ S Z 


— % + -dx-dy-dz + ■ dx- dy ■ dz+ ■ dx- dy ■ dz 


^=(-f + ^ + TS í + t 1 ) -dx-dy-dz 

dF z _ _ dP i d Ty Z dz xz i dx zz 

dx-dy-dz dz ' dy ' dx ' dz 


dF, = 

dv 





En total 3 components: 



dv 


* 4 Z — - |^+P*g+V • 
Jv <? z + V T z 


> 


y 


Quan tenim les 3 components una manera d’unificar-ho tot és posar un subíndex j i així ja deim que és un 
vector i és pot transformar en tres components. 
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Una densitat per la derivada de la velocitat respecta el temps: 



És equivalent a sòlids rígids 
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